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Abstract 



The question is studied whether weak solutions of linear partial integrodifferential equa- 
tions approach a constant spatial profile after rescaling, as time goes to infinity. The pos- 
sible limits and corresponding scaling functions are identified and are shown to actually 
occur. The limiting equations are fractional diffusion equations which are known to have 
self-similar fundamental solutions. For an important special case, is is shown that the 
i^h ' asymptotic profile is Gaussian and convergence holds in L 2 , that is, solutions behave like 

fundamental solutions of the heat equation to leading order. Systems of integrodifferential 
equations occurring in viscoelasticity are also discussed, and their solutions are shown to 
behave like fundamental solutions of a related Stokes system. The main assumption is 
that the integral kernel in the equation is regularly varying in the sense of Karamata. 
> 1 45K05; 35B40 
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CN ' 1. Introduction 



Consider the linear heat equation ut = Au in M. n , with fundamental 



solution U(x,t) = ^ ^n/2 e l^' 2 / 4 *. It is well-known and easy to see from 

^— > 

a 



the solution formula that as t —* oo, u(x, t) = U U(x,t) + o(t-™/ 2 ), where 
Uo = J Rn u(-, 0) is the initial mass of the solution, assumed to be finite. Thus, 
t n / 2 u(x\/i,t) — > UqU(x,1). Similarly, for solutions of the wave equation 
utt = u xx on R x [0, oo) with initial data u(-,0) = uo, u t (-,0) = 0, the 
well-known solution formula u(x, t) = ^ (uq(x + t) + uq(x — t)) implies that 
tu(xt,t) ~U \(5 _i + 5\) as t — > oo, in this case in the sense of distributions. 
For the case of the wave equation on M. n , the solution formulae that use 
spherical means imply that t n u(xt, t) — > Uow, where w is a distribution of 
dimension dependent order that is supported on the unit sphere in M. n . In 
the case of the heat equation, the solution depends (up to a multiplicative 
factor) asymptotically on the similarity variable £ = and convergence is 

uniform. In the case of the wave equation, the similarity variable is £ = -^r, 
with convergence in a space of distributions. 

In this paper, I investigate whether solutions of integrodifferential equa- 
tions 

ut(-,t) = ao Au(-,t) + / a(t — s)Au(-, s)ds 
Jo 

l 
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in R n have similar properties. Here R 9 ao > is a scalar, and a : [0, oo) — ► 
R is a scalar kernel. The special cases ao = 1, a = and oo = 0, a = 1 
correspond to the heat equation and the wave equation, respectively, and 
therefore are included. Thus the question is whether m(t)u(k(t)x,t) ~ Uoo 
as t — > oo, in a suitable sense. 

Since all equations in this class are of the form m + V • q = for some 
flux q, the ^-integral of solutions is formally preserved as t varies, 

u(-,t) = UQ = U( 



o • 

One therefore expects that m(t) = k(t) n , that is 

k(t) n u(xk(t),t) ~ UoWoo 

for a suitable function k(-) as t — > oo in a suitable distributional sense, for a 
limiting distribution Woq. Then the "trivial" behavior Woo = 5o can always 
be achieved by letting k grow to oo very rapidly. This trivial behavior must 
therefore be excluded. Also, solutions are expected to go to zero locally, so 
the trivial case Woo = is possible if k goes to oo too slowly and must also 
be excluded. With the right choice of k one hopes to obtain a nontrivial 
limit woo. 

In this paper, I show that for a large class of such integrodifferential 
equations there is a choice of k (unique up to an asymptotically constant 
factor) for which the limit Woo is indeed non-trivial. It turns out that the 
correct choice is 

k(t) = ^t^ciQ + J a(s)ds 

The limiting distributions Woo are also identified. They turn out to belong 
to a one parameter family, parametrized by (3 G (—1,1], with [3 = corre- 
sponding to the heat equation, [3 = 1 corresponding to the wave equation, 
and the cases of non-integer (3 corresponding to fundamental solutions of 
fractional diffusion equations. The main assumption is that the integrated 
kernel A(t) = ao + f a(s)ds should be regularly varying in the sense of 
Karamata ( 1J, and the index of variation (3 then determines the limiting 
distribution u^. As an aside, it should be noted that for the same Woo, there 
are many types of scaling functions k possible that are not asymptotically 
equivalent. It will be shown that all these possible limiting distributions are 
actually attained (in the sense of distribution, or in an important special 
case in L 2 ). 

The literature on self-similar asymptotics is huge, so I only mention the 
book |2j by G. Barenblatt. Fractional diffusion equations were discussed in 
[16] and [£j , with systematic studies carried out in El ■ Various 
physical models leading to fractional diffusion equations are discussed in 
I10| I18j . The main reference for integrodifferential equations of the type 
discussed here is the book jT^j by J. Priiss. The idea that regularly varying 
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integral kernels lead to asymptotically self-similar wave profiles for problems 
in viscoelasticity is exploited in ^3], for the case of the signalling problem. 
A related asymptotic concept is equipartition of energy, discussed for a class 
of exponential kernels corresponding to (3 = 1 in [3]. 

The plan of this paper is as follows. In the following section 2, two 
types of scaling (introduced at the end of this section) and the relations 
between them are discussed. All possible distributional limits and their 
corresponding scaling functions k are identified in section 3. Section 4 is 
devoted to giving sufficient conditions such that these distributional limits 
are actually attained. In section 5, the question of asymptotic self-similarity 
is studied in L 2 , leading to results about the time-asymptotic behavior of 
solutions that are sharp to leading order. The same question is taken up 
for three-dimensional linear homogeneous isotropic viscoelasticity in section 
6. Appendix A contains two important technical results for families for 
scalar integral equations, and two explicit examples for asymptotic behavior 
outside the theory developed in this paper are presented in appendix B. 

The notation (u, if) will be used for the result of applying a distribution 
u G V to a test function ip G V = C^°(M. n ). The pairing between test 
functions G C^°(R n x [0, oo)) and distributions U on R n x [0, oo) is denoted 
by {{U, <&)). In particular, if U(-,t)t>o is a family of distributions in T>' that 
is measurable and bounded with respect to the system of seminorms defining 
the usual topology on T>' , one can write 

/•oo 

(([/,*)) = / (U(;t),$(;t))dt. 
JO 

Convolution with respect to t G R is denoted by an asterisk, u * v(t) = 
J u(t— s)v(s) ds if u and v are both supported on the positive half axis. The 
Fourier transform of a function / G V is denoted by /, /(£) = L n e~ lx ^f{x) dx, 
and this is extended in the usual way to functions in L 1 or in L 2 or to dis- 
tributions. The Laplace transform of a function a : [0, oo) — > R is denoted 
by a(s) = / °° e~ st a(t)dt if defined, i.e. for s G C such that "Sis > a for 
some a G (— oo,+oo]. The usual notation for Lebesgue spaces L P (Q) and 
for the Sobolev spaces H s (M. n ) is employed, 1 < p < oo, — oo < s < oo. 
Vector-valued Lebesgue and Sobolev spaces are denoted in the usual way, 
e.g. LP([0,T],H 1 (R n )) or H S (R 3 , R 3 ). 

Let R + 9 t i — > u(-,t) be a measurable and locally bounded family of 
distributions on R n , and let k : [0, oo) — > R + be continuous and increasing 
to oo. In this paper, the scaled version of u is denoted by u^, defined by 

{u k (;t),<p) = (n(-,t),^(/c(t)- 1 -)) 

or, in case u(-,t) is a function for almost all t, 

u k (x,t) = k(t) n u(k(t)x,t) . 

Thus if u(-,t) G L 1 (R n ), then also u k (-,t) G L 1 (R n ), and the L 1 integral is 
unchanged. 
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There is an alternative scaling UT,k, denned by 

(u Tik (;t),<p) = (u(-,Tt),^(A;(r)- 1 -)) 

or if u(-,t) is a function for almost all t, 

u Tik (x,t) = k(T) n u{k{T)x,Tt) 

for x G R n and t > 0. The result now depends on k and the parameter 
T > 0. This scaling again preserves the L 1 -integral. 

2. Scaling 

Let 1 1 — > u(-,t) be a measurable and locally bounded (with respect to the 
usual system of seminorms) family of tempered distributions on 1", and let 
k : [0, oo) — > M + be continuous and increasing to oo. Assume that the scaled 
family Ufc(-,i) converges in T>' to some U G T>' as t —* oo, that is, 

(1) < Ujfc (. >t ) j¥ ,(.))_> (c/^v^ecniEr). 

It is necessary to relate this property to the behavior of the family UT,k as 
T — > oo. Using a tensor product argument, one easily sees that for all test 
functions $ G C^°{M. n x (0, oo)) 

(2) /°°K fe (- ; r), $(-, r))dr = (« fc (-, Tr), H^y; r))dr . 

Assuming that Q holds, the goal is to obtain a nontrivial limit in (j2J as 
T — > oo. Then a natural assumption is that 

< 3 > ^w =!,(T) 

exists for all r G M + , since in this case 

uniformly with all derivatives in x, boundedly in r, and thus 
(«r,fc(-,r),*(.,T))dr-> / (C/,$(p(r)-,r))dr. 



•X- 





A function k that is eventually positive and for which holds for some 
function p, for all r G (0, oo) is called regularly varying (|U). It is known 
that in this case p(r) = r a for some a G R which is called the index. An 
equivalent condition is 

for some b, p\ and all r in a neighborhood of r = 1. In this case, necessarily 
Pi(t) = Cr a for some C = pi(l) > 0, and one may choose b(t) = k(t). 

Since k is increasing b y assumption, a must be non-negative. If (jHJ) holds 
merely for r in a set of positive Lebesgue measure, it must already hold for 
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all positive r, and the limit is uniform on closed subintervals of (0, oo) (see 
PQ). Moreover, a can be recovered from the limit (which always exists) 

log k(t) 
a = lim — . 

t^co log t 

Returning to (J2J) and assuming now that k is regularly varying with index 
a > 0, one obtains 



(5) / (u T;k (;T),$(;T))dT -» / (U,<f> T (;T))dr 

Jo Jo 
where $ r (x, r) = <3?(xr a , r). If [/ is a locally integrable function, this means 

ur,fc(s,T) ~ r-" Q C/(xr- a ) . 

in a suitable sense (e.g. pointwise a.e. in (x,r)), as T — > oo. All this proves 
the first statement of the following proposition. 

Proposition 2.1. Let t u(-,t) be a measurable locally bounded family 
of tempered distributions on M. n . Let k : [0, oo) — > M+ be continuous and 
increasing to oo, and regularly varying with index a > 0. Define and 
UT,k as above. 

a) Lf U is a distribution on R n such that Uk(-,r) — > U, then UT.k(-,T) —* 
U(-, t) for a.e. r inT>' as T — > oo, where 

(6) (L7(-,r)^) = (C/,^(r a -)). 

Assume that U G -fP(R n ) /or some s£R and i/iai n-j^ converges to U, 
defined in®, locally uniformly in r as an i!P - valued function. Then also 

Uk(-,t) -> C7 

in # s (M n ), as t ^ oo. 

Suppose U G L r (M n ) /or some r € [l,oo] and Uk(-,t) — > U in L r (M. n ). 
Set w(x,t) = kityn-Uixkit)- 1 ), then 

\\u(;t)-w(;t)\\Lr=o(k(t) n ^ 1 -V). 

Moreover, if D m U E L r (M n ) and D rn u k (-,t) -> D m [/ in L r (R n ) /or some 
partial derivative D m of order m, then 

\\D m u(-,t) - D m w(;t)\\ L r = o(k(t) n ( r ~ 1 ~ r> ~ m ) . 



Proposition 2.2. Let k, I be regularly varying functions with index a > 
such that lim^oo = C E (0, oo), and let U E T>' . Suppose that 

u T,k(',t) ~~ > U(-,t), defined as in ©, locally uniformly in some H s with 
s G M. Then UTi(-,t) — * V(-,t) in H s , locally uniformly in t, where for 
ip G V 

(V(;T),i P ) = (U^(C' 1 T a .)). 

Thus if U is a function, then 

u T j(x,t) ~ C n T- na U(CxT~ a ) 
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in a suitable sense, and the estimates of part c) of the previous proposition 
hold. 

The easy proofs of the remaining parts of Proposition 2.1 and of Propo- 
sition 2.2 are left to the reader. 

If u(-,t) — > in some weak sense as t — ► oo, then it is possible to find k(-), 
going to oo, such that ux,k — > as a distribution on W l x (0, oo). Similarly, 
if u(-,t) G L 1 for all t and L„ u(-,i) = C is constant,then one expects that 
u T,k{'->t) C5q if &;(•) goes to oo sufficiently rapidly. The limiting cases 
U = and U(-,t) = C5$ should be excluded and will be called trivial. 
Non-trivial limiting distributions U should be neither zero nor supported on 
{0} x [0,oo) C W 1 x [0,oo). 

3. Identifying Asymptotic Limits 

In this section, I shall classify the types of limiting behavior that are 
possible for distributional solutions of the partial integrodifferential equation 

(7) ut(x,t) = aoAu(x,t) + / a(t — s)Au(x, s)ds 

Jo 

or more shortly Ut = aoAu + a * Au for x G M n , t > 0, with initial data 
u(-,0) = uq. It will be assumed throughout that ao > 0, the integral kernel 
a is bounded on any set [e, oo) and integrable on (0, 1), and no G L 1 (M n ). 
Let us begin by describing the limiting equations and their solutions. For 
a > 0, the Mittag-Leffter function E a is defined as 

oo ^ 

see 0. This is an entire function for any a > 0. Special cases include 
E\(z) = e z , E2{z 2 ) = cosh(z), and Eu 2 (z) = e z erfc(-z), where erfc is 
the complementary error function. For < a < 2, a ^ 1, there is an 
asymptotic expansion 

N -n 

£rra + 0(|zr "' ) 

71=1 ' 

as z — > oo in a sector about the negative real axis, where the reciprocal of 
the r - function is extended as zero at the poles of T. In particular, E a is 
bounded on the negative real axis for all a < 2; see jS] for details and other 
properties. 

For a > 0, A G R the function u(t) = E a (—Xt a ) is the solution of the 
scalar integral equation 

(9) «(*) + rTT [\t-s) a - 1 u(s)ds = l 

"(a) Jo 

as a direct calculation shows. 
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Let < a < 2. Define w a (-,t) as the tempered distribution on W 1 whose 
spatial Fourier transform is 

(10) tia(H,t) = E a (-\t\H a ) (£ G M n , t > 0) 

for t > 0. If a < 2, the asymptotic behavior of -E a implies that w a (-,t) G 
.£P(IR n ) for s < 2 — ra/2. In particular, if a < 2, with the exception a = 1, 
then w a (-,t) G L 2 (R n ) if and only if n < 3andw Q (-,i) G L 00 ^™) if and only 
if n= 1. For a = 2, one has u> 2 (£,*) = cos(|^|t) and thus w 2 (-,t) G iT s (K n ) 
for all s < — n/2. If a = 1, then w)i(£,i) = e - '^' *, and w± is the well-known 
fundamental solution of the heat equation. Also, w a (-,t) —* <5o as i — ► 
for all a, in the sense of distributions. 

The distribution w a solves the integrodifferential equation 

(11) w(-,t) - A (J r(a)^ 1 (t - s^wi-, s) ds^j = 5 . 

This follows immediately from @ and (jlUj) . If a = 1, this is the heat 
equation, and if a = 2, this is the wave equation. For 1 < a < 2, equation 

(12) can be differentiated formally, resulting in the fractional heat equation 

(12) w t (-,t) = A (J* T(a - l)~\t - s) a - 2 w(; s) ds\ . 

These equations have been studied in pi 171 IUI fTTl [T2l [Tb], 

Returning to solutions of (|7|). let us assume that i i— > u(-,t) is a func- 
tion with values in the set of tempered distributions such that for all test 
functions $ G Co°(M. n x [0, oo)) the equation holds 



/ «,s),$ t (-,s)+a Q A$(-,s)+ / 

j0 



a(t - s)${-,t)dt) ds 

+ M,0),$(-,0)> =0. 
With = ao + J one can write equivalently 

(13) «(-, t) = u(-, 0) + A A(t - s)u(-, s) dsj 

in the sense of distributions, or with &±(x, s) = A(t — s)<J?(cc, t) dt 

/■OO /'OO 

(14) / a) -A*!(., a)) ds = («(-, 0), / $(-,t)eft>. 
jo io 

Let us consider solutions of (J7J) with the scaling 

wr, k {x,T) = k{T) n u{k{T),TT) 

introduced earlier, where k is left unspecified for now. Set K = k(T), then 
v = u T,k is seen to be a distributional solution of the problem 

T 

(15) V = U 0jK + -j£2 A T * Av 
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Tt 

where uq^k = K u uq{Kx) and Ax(t) = A(Tt) = a$ + f a(s)ds. Also, the 
spatial Fourier transform ut k solves 

(16) u T Att) + m 2 ^A T *u T ^,t) = M^) 

in the sense of distributions. If one wishes to obtain a limiting equation of 
a similar form, one is led to assume that there exist functions p, A^ such 

lim ^=A O0 (t). 

T^oo p(T) 

Let us also assume that A is eventually positive (not necessarily bounded 
away from zeros). As explained in the previous section, this implies that A is 
regularly varying and that one may choose p(T) = cA(T) = cAt(1), Ax>(*) = 
c-V for some j3 £ M and any constant c > 0. If the kernel A^ is to be 
integrable at t = 0, then one should require > — 1. Since vl' = a was 
assumed to be bounded on (l,oo), necessarily (3 < 1. This motivates the 
main assumption in the following result. 

Theorem 3.1. Let u be a solution of (J7J) in the sense described above, and 
let A be eventually positive and regularly varying with index j3 S (—1,1]. 
Assume that there exist a non- decreasing function ko : [0, oo) — ► [0, oo) with 
ko(oo) = oo, a sequence T n — > oo and a non-trivial limiting distribution Woo 
on K n x [0, oo), Uoo(-,t) £ H s (M. n ) for a.e. t for some fixed s > —oo such 
that 

u T n ,k —> u oo 

as n — > oo, a.e. boundedly in H s (R. n ). Then one may choose k(t) = 
\J tA(t)T(l + p), and with this choice 

UT,k -> U W 1+l3 = Uoo 

where Uq = J Kn uq(x) dx and w\ + p is the distributional solution of the inte- 
grodifferential equation 

w(-, t) - A (J* r(l + P)~\t - 8 fw(; s) ds^j = 5 

defined in ()1U[) . 

Proof. It should be noted that A is regularly varying with index P iff k is 
regularly varying with index The assumptions for A imply that 



A(T)T{1 + P) 00 w T(l + P) 
as T — > oo. Let $ G Co°(R n , x(0, oo)) be a test function such that 

((A^oo.Uoo)) / 

/■oo 

((««>,*» / / $(0,t)dxdt 

JO 
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where <3?oo (x, s) = j4oo(* — s)$(x, t) dt. This is possible since the limiting 
distribution is non-trivial (not identically equal to zero, not supported on 
a subset of W 1 x [0, oo) and in some fixed H s for a.e. t. Set <& n (x, s) = 
It° ^^"r(i+/3) ^( x ^) dt. Then <3? n — > together with all derivatives. 

Let us first show that ko may be replaced with k, i.e. ux n! k — * '"oo- 
Indeed, since 

«*,«r n ,*,»= /°°(*(-, tUoMT n )) dt + r(1 + jf)^ Tn) «A* n , n TnM )) 

JO Ko{l n ) 

and all terms except the fraction have limits as n — > oo, it follows that 

L = Hm r(l + (5)T n A{T n ) 

n^co k (T n ) 2 

exists. Of course, L > 0, and after replacing ko with one may assume 

without loss of generality that L = 1. One can therefore replace ko with 
k(t) = y / r(l + /?)L4(t) and obtain that UT n ,jfc ~~ ► ""oo- 

Observe next that ()16|) now takes the form 

I6I 2 £ 
(17) u Tife (C, t) + A(r) r 1 + j9) 4r * *r,*& «) = fio(^) 



"T,fc(£,*) = z T(\t\ 2 ,t)u [ — 



for all £. Thus one can write 
where Zt(p, ■) solves the equation 

(18) ^') + ^ iT " T(ft ' ) = L 
Note that zx(p,t) = z(X,Tt), where z(A, •) solves 

(19) z(\,-) + \A* z(\,-) = 1 

WWII A — rj 4(T)r(i+/3) • 

By Lemma A. 2, as T — ► oo, Zt(p, t) converges to E 1+ p(-pt 1+f5 ), locally 
uniformly in p and t > 0, and thus «T,fc converges pointwise in £, locally 
uniformly in t, to v{^,t) = UoEi + p(— \i\ 2 t l+l3 ) = £^o^i+/3(£) i), that is, a 
solution of the limiting equation 

P 



But for the subsequence T n , the limit is Uqq. Therefore, convergence holds 
along the full sequence T — ► oo, and the limit is UoWi + p. By Parseval's 
identity, the theorem follows. □ 
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To sum up, the possibilities for limiting behavior identified in this result 
are the following: 

1. Behavior like the fundamental solution of the wave equation {[3 = 1), 
expected if e.g. a(-) ~ c > or a(t) ~ (logi)" 1 as t — > oo for some real 
number m. In this case, k(t) ~ y/ct or k(t) ~ t(logt) m / 2 . 

2. Behavior like the fundamental solution of the heat equation (/? = 0), 
expected if e.g. a(-) is integrable and either a > or a(s)ds + ao > 0, 
but also if e.g. a(t) ~ i -1 . If a is integrable, then k(t) ~ V^4i, where 
A = ao + / °° a(s)ds, while if e.g. a(i) ~ i -1 , then ~ y/t logi. 

3. Behavior like the fundamental solution of a fractional integrodiffer- 
ential equation of order 1 + (3. If < /3 < 1, this is expected if e.g. 
a(t) ~ t 13 ^ 1 . If — 1 < ft < 0, this may occur if e.g. a is negative and inte- 
grable, J °° a(s)ds = —ao, and a(t) ~ — t@~ l . In each case, k(t) ~ t( 1 +/ 3 )/ 2 . 
Of course, the behavior of k may be modified by additional logarithmic 
factors also in this case. 

4. Asymptotic Distributional Limits 

The purpose of this section is to demonstrate that under mild assump- 
tions, solutions of (7) do converge to limiting solutions under the scaling 
u ~» iPrk- I shall start by stating a general existence result for strong solu- 
tions that is essentially well-known for the Hilbert space case; see |15j . 

Proposition 4.1. Let uq G H s (M. n ) for some s£l and assume that 

lim(a + ytaiiuj + h)) > 

/or a// a; G M. T/ien t/iere exists a unique function u £ C([0, oo), H s (M. n )) 
that solves Q in t/te sense of distributions and for which u(-,0) = uq. The 
Fourier transform u is given by 

(20) u(tt) = z(\a\ 2 ,t)u (Z) 

where z(X, •) is the solution of (fl9|) . 

The condition for the kernel a is equivalent to the requirement that its 
cosine transform is bounded below by — ao in the sense of measures. Equiv- 
alently, the measure ao<5o + a(t)dt is required to be positive definite (see [B]). 
We do not require that uq G L 1 and cannot assert that u(-, t) G L 1 for t > 0. 

Proof. Consider the integral equations (|19|) for A > 0. By Lemma A.l, the 
estimate |z(A, t)\ < 1 holds for all A and t. Then define u(-,t) as in l|20|) . 
Using Parseval's identity and the bound for z(X, •), one obtains \\u(-, t)\\H s < 
||^o||h s for all t. Also, u(^,t) — ► «o(0 pointwise a.e., and by construction 

^(e, t) + |£| 2 (a n(e, t) + a * t)) = 

for almost all ^. Since 2f(A, •) is continuous, locally uniformly in A, and 
uniformly bounded, u(£, •) is also continuous with values in H s (W l ). Thus u 
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is a distributional solution of l[7]). Uniqueness follows by taking the Fourier 
transform of a solution in this class and recognizing that it must have the 
form ((201). D 

Such solutions converge to the limiting solutions identified in the previous 
section under the scaling introduced there, if the primitive A of a is regularly 
varying. As explained above, these are natural conditions. From now on uq 
will always be assumed to be integrable. 

Theorem 4.2. Assume that 

1. the kernel satisfies linnjjo( a o + ^o,{iu) + h)) > for all u £ M. 

2. the primitive A(t) = ao + J* * a(s)ds is eventually positive and regularly 
varying with index (3 £ (—1, 1]. 

Set k(t) = y/tA(t)T(l + pi), then for the solution u of (|7|) found in the 
previous proposition 

(21) «r,k(-,t) -» U oUoo (;t) 

in H s (R n ), if s < -n/2. Here = E 1+(3 (-\£\ 2 t 1+l3 ) and U = 

J^n u o(x)dx. The convergence is uniform on any interval [c,d] C (0, oo). 

A few remarks can serve to put the result in perspective. First, if P = 1, 
then Uoo(C,t) = cos(|C|i), and thus Uoo(-,t) ^ H~ n / 2 (R n ). Thus for a result 
that covers the entire range P £ (—1,1], one cannot expect convergence in 
better spaces than H s (W n ) with s < —n/2. Also, conditions 1 and 2 in the 
above result are independent. For example, the kernel a(t) = cos(t) with 
A(t) = sin(t) and a = \ (S{ + 5-i) in the sense of measures on iM. satisfies 
condition 1, but is not regularly varying. The kernel a(t) = 1 — e~ l has the 
antiderivative A{t) = ao + t — 1 + e~ l which is regularly varying with P = 1, 
but since ^fcd(iuj) = — 1 ^ for uj ^ 0, condition 1 is not satisfied if ao < 1. 
Finally, it should be recalled that L l {R n ) C H s (R n ) for s < -n/2, but not 
for larger s. 

Proof. Let u be the distributional solution in C([0, oo), H s (U. n )) constructed 
in proposition 4.1. Then ux,k satisfies 

|£| 2 / ( \ 

UT, k {d,t) + A{T)ni + p) A T * u T , k {i, t) =*>[ m ) 

and therefore with k(T) = ^TA(T)T(1 + P), 

where z(X, •) solves (|19p and therefore v\(t) = z ^ TA(T)h{i+i3) > ^) s °l ves 
(|1T7|). By Lemma A.2, as T -> oo, 
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for all A > 0. Consequently, 

u T ^,t) ^ u (0)E 1+(3 (-\^\ 2 t 1+ ^ 

pointwise for all £, locally uniformly in t. 

To prove that convergence holds in H s (M n ) for s < — n/2, one invokes 
again Lemma A.l to deduce that |#(A, t)\ < 1 and therefore 

\u T ,k(^t)\ < l«o (my 1 - c 

for all £ and T. Then for t > 

IK fc (-,t)-£W-,t)ll!p = / (l + \C\y\u T .k(^t)-U E 1+p (-\^t)\ 2 . 

Since 2s < — n, Lebesgue's dominated convergence theorem implies the con- 
clusion. 

□ 



5. Asymptotic Limits in L 2 

In this section, equation (J7J) will be considered under the scaling u «T,fc> 
with the goal of proving that limiting solutions are attained in L^ c (0, oo; L 2 (W 
or more generally Lf£ c (0, oo; H s (M. n )) with s > 0, provided the initial data 
are in 1? or H s . The limiting solutions were identified in section 2 and are 
unbounded in any L r (M n ), r > 1 as t —* 0. Thus one cannot expect uniform 
convergence up to t = 0. A convergence result in L 2 or in a better space 
allows one to obtain the exact asymptotic behavior of solutions of (7) in this 
space to leading order, by Proposition 2.1(c). 

I only have a result for the case where A is regularly varying with index 
(3 = 0, i.e. < a + / °° a(t)dt < oo. The result is independent of the space 
dimension. In this case the limiting equation is the heat equation. Note 
that for (3 ^ 0, the limiting distributional solution is in L 2 (M n ) for t > if 
and only if n < 3. Thus a result that holds for all spatial dimensions cannot 
be expected if /3 ^ 0. 



Theorem 5.1. Assume that 

1. the initial data satisfy u G L 1 (M n ) n i/ s (M n ) /or some s > 

/or some a > 0, a can be extended to the half plane {s G C|5Rs > —a} 
and either ao > and ao + a(iLj)) > 0, or 5fta(— a + iw) > for all ui. 
Set k(t) = ^tA{t), then for the solution u of (J7J) found in Proposition 3.1 

(22) «r,fc(-,t) -> I7o«oo(-,*) 

m iiP(R n ), w/iere f7o = J K „ ^(a^cte. Here Uoo(£,t) = e~l^ 2 * ; 

i/te fundamental solution of the heat equation. The convergence is uniform 

on any compact subinterval of (0, oo). 
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Proof. Note that assumption 2 implies that < ao + J °° a(t)dt < oo. Thus 
A is regularly varying with index f3 = 0. From the proof of Theorem 4.2 one 
sees that 



(23) u Tk ^t) = zU^,Tt)u 



fc(T) 



z(A,Tt)u 



fc(T) 



with A 



1 1|2 

t ^/ T n ■ By Lemma A.l, there are therefore estimates, valid for all 
sufficiently large T, 

<2e - C |?l 2 * 



if |£| 2 < dT and 



\TA{TY 



< 2e 



if |£| 2 > dT. Here c 
as in Lemma A.l. Now consider 

\\uT,k(-,t) - ?7(yUoo(-,*)l|H s 
(l + |£| 2s )|z(A,Tt)n 



< 



+ 



+ 



d = are positive constants, and e, L are 



k(T) 



[ (i + \a\ 2s )\ z (\,Tt)u 

J\£\<R 

(i + \a\ 2s )\ z (\,Tt)u 



iei>fl 



fc(T) 
fc(T) 

1/2 



\ 1/2 

I7 e-lfl 2 *| 2 de) 
- U Q e-\t\ 2t \ 2 dZ 

\ 1/2 



1/2 



.'\a\>R 

for arbitrary i?, where as before A 



(1 + \i\ 2s ) U 2 e- 2 ^ 



TAjTj- Then for given 5 > and [a, b] C 
(0, oo), choose R so large that the last integral is less than 5, uniformly in 
t £ [a, b], and then choose T large enough such that the first integral is also 



TA(T) 



,Ti) 



bounded by S, uniformly in t. This is possible because z 

e ~l€l 2 ' as T — > oo by Lemma A. 2, locally uniformly in £, and because uq is 
continuous. Then if dT > R 2 , the second integral in the last expression can 
be estimated by 

\ 1/2 



■■■<([ (i + iei 2s ) 

\J\t\>R 



e- 2cdTt \u 



k(T) 



2 ^ 



< e- cdTt k(T) n / 2+s \\uo\\ H s 



which is also smaller than S, if T is chosen sufficiently large, locally uniformly 
in t. The proof of this theorem is therefore complete. □ 



Using Proposition 2.2, one obtains 
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Corollary 5.2. Under the assumptions of Theorem 5.1, the solution u of 
satisfies 

\\ U (;t)-w(;t)\\ H s = (W 4 ) 

where w(x,t) = Uq (4:A OQ Trt)~ n ^ 2 exp (— \x\ 2 / (4A oa t)^ is the solution of the 
heat equation 

w t = A 00 Aw, w(-,0) = U 5 
and = an + J °° a(s) ds. 

6. Linear Viscoelasticity 

Consider a viscoelastic material with mass density p = 1 occupying all 
M 3 , and denote the displacement of a material point at position x and time 
t by u(x,t) and the velocity at this point by v(x,t) = dtu(x,t). Let us 
assume that the material is at rest for t < and prescribe an initial velocity 
field v(-,0) = vq. The Boltzmann model for linear isotropic homogeneous 
viscoelasticity (^S]) leads to the equations of motion 

(24) v t = a Av + — ^VV • v + a * Av H * VV • v . 

3 3 

Here an > 0, 6o £ M, and a, 6 are suitable scalar-valued functions that 
describe the stress response of the material under shear and compression, 
respectively. 

The reader should recall the well-known decomposition into divergence 
free and gradient components, as follows. For u G .£P(IR 3 , M 3 ), let Pu, Qu G 
.fP be defined by 



Then for an ff s -valued solution v of (|24|). one obtains that p = Pv and 
q = Qv satisfy the equations 

(25) Pt = (3 Ap + (3*Ap 

(26) q t = a Aq + a * Aq 

with /? = 4a °+ 2b ° and /3(t) = 4aW + 2fe(t) and with initial data p(.,0) = p = 
Pvq and g(-, 0) = qo = Qvq = vq — po- In addition, Vxp = and V • q = 
in the sense of distributions. Thus p and g satisfy scalar integrodifferential 
equations, and as in Proposition 4.1 one obtains an existence result together 
with a representation formula for the solution, given next. 



Proposition 6.1. Let vq G if s (R , R 3 ) for some s£l and assume that 

(27) lim(a + ~Std(iuj + h)) > 

hio 

(28) lim(A) + + h)) > 
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for allu G R. Then there exists a unique function v = p+q G C([0, oo), H S (M. 3 , 
that solves (|2l|) in the sense of distributions and for which v(-,0) = vq. The 
Fourier transforms p, q are given by 

(29) p{£,t) = ^i(|C| 2 ,t)po(0 

(30) g(£,t) = z(|£| 2 ,i)<7o(£) 

where z(X, •) is i/je solution of (|19|) and zi solves (|19[) wzi/i ao, a(-) replaced 

Let us now consider the case where a, /3 G L , corresponding to a vis- 
coelastic material with vanishing elastic equilibrium response, i.e. a liquid. 
In this case, the limiting behavior is expected to resemble the fundamental 
solution of the compressible Stokes system 

(31) w t = AAw + {B- A)VV ■ w 

where A = ao + Jq°° a(s) ds and B = (3q + J °° /3(s) ds. The fundamental 
solution is known to be the matrix valued function W(x,t) = U(x,Bt) + 
V(x, At), where 

with E denoting the identity matrix. In real terms, U and V can be ex- 
pressed in terms of error functions (Kummer functions, confluent hypergeo- 
metric functions), e.g. 

f«M = M-(^/(|L)). 

A recent derivation of these fundamental solutions in a more general situa- 
tion may be found in Note that these functions are not integrable with 
respect to x, since their Fourier transforms are not continuous at £ = 0; in- 
deed they behave like 0(|x| -3 as \x\ — ► oo for fixed t > 0, due to well-known 
asymptotic results for the Kummer function. Using the arguments that led 
to the proof of Theorem 5.1, one can now describe the asymptotic behavior 
of solutions of (|2l|) in terms of U and V. For this purpose, let us assume 
the following: 

• For some s > 0, v G L 1 (M 3 ,1R 3 ) n H S (R 3 , R 3 ). 

• For some a > 0, the Laplace transforms d, (3 can be extended to the 
half plane {z \ > — a} 

• Either ao > and ao + a{iuj)) > 0, or 3f?a(— a + iuj) > for all u. 

• The same assumption for (3q and /?(•). 

As in the previous section, one can then use the representation formulae 
(|29JI and (|30|) together with the results of Appendix A to prove the following 
result. 
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Theorem 6.2. Under these assumptions, the solution v of (|24jl satisfies 

\\ V (;t) - V Q T U(; Bt) - V Q T V(; At)\\ H s = (W 4> 



where Vq = vq(x) dx G M 3 , U and V are the components of the funda- 
mental solution of the compressible Stokes system (j31|) . and 

A = a + J a(s)ds, B = -A + - ^6 + J b(s)dsj . 

The result shows that to leading order for large t, the solution v(-,t) of 
the Boltzmann system (|2"4"|) behaves like the solution of the Stokes system 
(|31|) with distributional initial data w(-,0) = Vq5q. 

Appendix A. Scalar Integral Equations 

In the following, let us assume that z : [0, oo) — > M is a solution of the 
scalar integrodifferential equation 

(32) z'(t) + A (a z(t) +a* z{t)) = 0, z(0) = 1 

where a > 0, A > 0, and a £ Lj oc (0, oo;R), a G L°°(l,oo;R). Let a be 
the Laplace transform of a, defined for > 0. Recall that by Parseval's 
identity, ao + $ta(s) > c for all s in the right half plane, for some non-positive 
constant c, if and only if 

/ u(t) (a u(t) + a * u(t)) dt > c \u(t)\ 2 dt 
Jo Jo 

for all real- valued square integrable functions u and all T; see (Hj. 

Lemma A.l. 1. If limMo(«o + $ta(iu + h)) > /or all u € R t/ien 

\z(t)\ < 1 

/or a// 1 and all A > 0. 

Assume that for some a > 0, a can be extended to {s G C|5Rs > —a} and 
that either oq > and inf w6 ig(ao + a(iu)) > 0, or that 5Ra(— a + iw) > for 
all u). Then there is a constant e > such that for all t > 

(33) \z(t)\ <2e- €minix ' 1)t . 
Proof. Let us consider the more general equation 

(34) z'(t) + \(a + a*z(t)) = f(t) 

where / G -kj" oc (0, oo; R). It will be shown that 
a) under the assumptions of part 1, for A = 1, 



(35) 



K*)I<1+ f\f{s)\ds 
Jo 



for all t, 

b) under the first set of assumptions in part 2, with / = 0, 
(36) \z(t)\ < e~ 5t 
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for all t, 

c) under the second set of assumptions and with / = 0, 

(37) \z{t)\ < 2e~ st 

for all t. Here 5 = A min{e, 1}, and e > depends only on ao and a. Together 
these assertions imply the lemma. 

To prove part a), multiply ()32|) with z(t) and integrate over [0, T], result- 
ing in the identity 

\\z{T)\ 2 + f z(t) (a z(t) + a * z{t)) dt = \+ [ z(t)f(t)dt . 
1 Jo 1 Jo 

Since ao + d(ico) > 0, the integral on the left is non-negative, and the in- 
equality ||z(T)| 2 < \ + J |z(i)||/(i)|di follows for all T. Bihari's theorem 
now implies (|H5j) . Evidently this estimate is independent of A > 0. 

For the proof of b), set zg(t) = e St z(t) and as(t) = e St a(t), where S > 
will be fixed later. Then z$ satisfies 

(38) z' 5 (t) + A ((a - oA" 1 )^) + a s * z 5 (t)) = 

and zs(0) = 1. Note that ds(s) = d(s — 5) whenever 5 < a. Since 3f?d(s) is 
bounded on any vertical line = (5 with (3 > —a, harmonic to the right of 
any such line, and bounded away from near $ts = —a one can find e > 
such that ao + 3fta(— e + iuj) > e for all u> £ M. Then also ao + 3fta(s) > e 
whenever 5Rs > — e. 

Let now A < 1. Set S = Ae, then obviously 

ao — SX^ 1 + $las(iuj) = ao — e + $ld(—5 + iuj) > . 

If A > 1, one sets 5 = e and obtains 

ao — 5A _1 + $lds(iuj) > ao — e + 3f?a(— e + iu) > . 

Part a), applied to (j3*H|) . implies the desired estimate in both cases. 
To prove c), note that 

(39) MaUoj) > Cl „ 

LO Z + OL 

for all u>, for some c\ > 0, since 9ftd is harmonic and positive for ?R.s > —a. 
Also, argd(z) = for ?R.z = and arga(— a + iuj) > — f for u > by 
assumption. Therefore by the maximum principle for harmonic functions, 
argd(z) > arg(z + a) for all z with Qz > 0, ?R.z > —a. This implies that 

Qd(iuj) > —auj$id(iu)) 

for all uj > 0. Now let b(t) = de~ at with 

where c\ is as in ()39|) . Thus 5(0) = d, b' = —ab, and b(s) = Consider 
the kernel 

ai(t) = a(t) + Xa * b(t) - ab(t) 
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with Laplace transform 

~ f \ , A Xd \ ad 

aUs) = a(s) 1 H 

w w \ s + a) s + a 

for < A < 1. Then d\ is analytic for 5is > —a, and for s = iu>, to > one 
has 

m~ /■ \ <x\~ i ■ \ ( i ^da \ .. . Xdto a 2 d 

;Rai{iuj) = }Ra{iuj) 1 H 7, ^ + zsa{iu)- 



uj 2 + a 2 J to 2 + a 2 to 2 + a 2 



2(w 2 + a 2 ) 2 ; a; 2 + a 2 u 2 + a 2 

> 

uniformly in A G [0,1], by the choice of d. In addition, Kai(s) is bounded 
uniformly in A and s on the strip —a/2 < !Rs. One can therefore find a 
positive 7 < a such that whenever 3fts > —7, then | + 3?ai(s) > 0. 

Now the estimate (|37[) can be proved for small A, say A < min{l, |^}. 
Forming the convolution of (|32j) with Xb and adding the result to (|32|) . one 
obtains the equation 

z'(t) + Xb * + A (a + A6 * a) * 2 = 

or equivalently 

z'(t) + A (dz{t) + ai * z(t)) = Xb(t) = Xde~ at . 

where a\ is as above (depending also on A). Set 5 = 4p and as before 
zs(t) = e St z(t), a^sif) — e St ai(t). The resulting equation for z$ is 

z' s (t) + A (p s (t) + ai, s * z s (t) \ = Xde^\ z s (0) = 1 . 

Since $lai t $(iu}) = $ldi(—5 + iuS) > — |, part a) of this proof implies that 

\z s (t)\ < 1 + / Ade^+ Ad '/ 2 = 1 + — ^— < 2 

where the last inequality follows from Xd < I7. This implies 

\z(t)\ < 2e~ Xet 

with e = f , whenever A < min{|^,l}. Reducing e if necessary, inequality 
(|H7|) is proved for A < 1. 

The proof for A > 1 is similar. One considers the kernel function 

a 2 (t) = a(t) + a * b(t) - A _1 a6(i) 

with Laplace transform 

d \ X~ 1 ad 



a 2 {s) = d(s) 1 + 

s + a J s + a 
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Then by the same argument, | + 5Ra2(s) > whenever "Sis > —re, for some 
k > that does not depend on A > 1. Then form the convolution of (|32|) 
with b and add the result to (|32|) . This yields the equation 

z'(t) + dz(t) + A (a 2 * *(*)) = b(t) = de~ at . 

with a 2 {t) = a(t) + a * b(t) — A~ 1 a6(t). Set 5 = min{re, |,a — d} and 
z s = e st , a 2 ,s = e St a 2 (t). Then 

z' s (t) + (d- 5)zs(t) + A (a 2 ,s * z s (t)) = de^*, z 5 (0) = 1 . 

Since ?R.d 2 s(iuj) = $ld 2 (— 5 + iuj) > — |, part a) again implies that 

r-oo J 

\z 5 (t)\<l+ de^ * = 1 + — — <2, 
Jo « - o 

where the last inequality follows from 5 < a — d. Therefore, 

\z(t)\ < 2e~ 5t < 2e~ et 

if e as chosen earlier or possibly lowered, whenever A > 1. The proof is now 
complete. □ 



Lemma A. 2. Let (A n ) n >i be a sequence in L 1 (0,To;M) such that \\A n — 
AxdIIl 1 — * 0- F° r P > 0; let w n{p, •) be the solution of 

w n (p, t) + pA n * w n (p, t) = 1 

for0<t < T . Then 

W n (p, •) ->■ Woo(p, •) 

uniformly in t £ [0,Tq], locally uniformly in p, where 

Woo(p,t) + pAoo * Woo(p,t) = 1 . 

In particular, asume that A is a regularly varying kernel with index (3 £ 
(—1,1] and eventually positive, and set Ax(t) = A(Tt). Then the solutions 
vx(-) of 

(40) v ^ + MfWW) AT * v ^ = 1 

converge uniformly in t and locally uniformly in A to Ei + g(— \t 1+ @). 

Proof. The first assertion follows from a standard argument for Neumann 
series. Since ^xyr(i+/3) ~~ ^ r(i+/3) ' P orn t w ise in t and also in L 1 (0,To), the 
second assertion follows as well. □ 
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Appendix B. Two Examples 



Here are two explicit examples of integrodifferential equations of the form 
for which the assumptions in the main results are not satisfied and the 
conclusions fail as well. 

First consider (J7J) with the kernel a(t) = cos(i) and ao = 0. Thus a(s) = 
pq-j-, and therefore this kernel is positive definite; ?R.a(s) > for all s in the 
right half plane. Since A(t) = sin(t), the kernel is not regularly varying. 
Taking the Laplace transform with respect to t and the Fourier transform 
with respect to x, one obtains that the Fourier-Laplace transform solution 
u satisfies 



After solving for u and inverting the Laplace transform, one obtains 



Thus u(x,t) = u\(x) + U2(x,t), where ui — Aui = uq and U2 solves the 
Klein-Gordon equation U2,tt + U2 = Aii2 with initial data U2(-,0) = uq — 
ui, U2 ) t(-, 0) = 0. Locally in x, u(-,t) — > u\ as t — > oo, since the contributions 
from U2 are radiated off to infinity. There is a nontrivial time-asymptotic 
limit (attained e.g. pointwise a.e. for sufficiently smooth initial data) that 
depends on the initial data. 

As a second example, consider Q with the kernel a(t) = — e _< and ao = 
1. Thus ao + a(s) = and this kernel is also positive definite. Here 
A(t) = e , and the kernel A can be viewed as regularly varying with index 
(3 = — oo. As before, taking the Laplace transform with respect to t and the 
Fourier transform with respect to x, one obtains that the Fourier-Laplace 
transform solution u satisfies 



The equation can be solved for u and the Laplace transform can be inverted, 
and the result is 



In this case therefore u(x, t) = u\{x)+U2(x, t), where as before u%— Au\ = uq 
and ii2 now solves the diffusion equation u-it + u 2 = Au2 with initial data 
^2( - ,0) = no — u\. Again, there a nontrivial time-asymptoptic limit that 
depends on the initial data, namely u(-,t) — u\ = 0(t~ n / 2 e~ t ) as t — > oo, 
uniformly in x. 




u = Uq . 





u = Uq . 
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